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We compute one-loop correlation functions for the fluctuations of an interface using a field theory 
model. We obtain them from Feynman diagrams drawn with a propagator which is the inverse 
of the Hamiltonian of a Poschl- Teller problem. We derive an expression for the propagator in 
terms of elementary functions, show that it corresponds to the usual spectral sum, and use it 
to calculate quantities such as the surface tension and interface profile in two and three spatial 
dimensions. The three-dimensional quantities are rederived in a simple, unified manner, whereas 
those in two dimensions extend the existing literature, and are applicable to thin films. In addition, 
we compute the one- loop self-energy, which may be extracted from experiment, or from Monte Carlo 
, simulations. Our results may be applied in various scenarios, which include fiuctuations around 

' topological defects in cosmology, supersymmetric domain walls, Z(N) bubbles in QCD, domain 

' walls in magnetic systems, interfaces separating Bose-Einstein condensates, and interfaces in binary 

, liquid mixtures. 
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I. INTRODUCTION 

Many natural systems exhibit interfaces that separate regions of different physical characteristics. An interface in 
, a binary liquid mixture separates its two components, a domain wall in a magnetic system separates its magnetic 
' phases. Apart from those traditional examples, several systems of current interest can be viewed as different guises 
. of that same physical situation. They include topological defects in cosmology }Tj, supersymmetric domain walls 1^, 
■ Z(N) interfaces in thermal SU(N) gauge theories^, or different types of Bose-Einstein condensates 4]. 
, Recently, studies of quantum and statistical fluctuations around interfaces or domain walls have been the object of 
• renewed attention. Such studies have concentrated on one-loop calculations which, technically, amount to computing 
fluctuation determinants around the interface (domain wall) background. They include computations around kink 
'" ^ ^ backgrounds in scalar theories in various dimensions 'sj, as well as in supersymmetric models (^. They use different 
I ' methods, exploiting connections with special properties of determinants of differential operators^, with scattering 
Qh, data@, and with the spectrum of the operators^- 

^ ' The calculations mentioned in the previous paragraph are restricted to vacuum bubbles, as will become clear in 
'T^ \ the sequel. In the present article, we will go beyond by computing correlations that involve one and two-point 
^ functions. Determinants, vacuum bubbles, and correlations will all be obtained from a semiclassical propagator that 
describes how the fluctuations of an interface (domain wall) evolve. Our method relies on a systematic semiclassical 
^ expansion around the given background, and has the semiclassical propagator as its essential ingredient. It not only 
' serves as an alternative to the methods used in particle physics inspired applications HSH, but also provides a 
unified framework that extends those results to the computation of correlations, allows for connections with statistical 
mechanical systems (for which one-point functions have been obtained by other methods), and introduces calculations 
of two-point functions that lead to novel results. 

We shall profit from the connections with statistical mechanical systems to present our method. Indeed, it is well- 
known0, Q that interface fluctuations may be described by a scalar field theory model with a double-well potential. 
The model admits a classical solution, a kink profile depending on only one spatial (longitudinal) coordinate, which 
is associated with the mean-field configuration of the interface. 

Fluctuations of the interface have been taken into account in 4 — e dimensions, via renormalization group methods0, 
0, as well as directly in three-dimensions |0, a result of their inclusion, calculations of the modified surface 

tensionj3, 13 1 and of the modified interface profilell 113 were successfully carried out: the surface 

tension was computed up to two-loop orderp^. leading to the prediction of universal ratios for an interface in the 
three-dimensional Ising model (which belongs to the same universality class of the scalar field model); the interface 
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profile was computed up to one- loop order |9|,lll|, leading to a detailed comparison with experimental data[l5j for the 
reflectivity of binary liquid mixtures near a phase transition. 

The three-dimensional calculations listed above can all be viewed as the result of computing the Feynman diagrams 
of a semiclassical expansion around the mean- field interface background T^. The diagrams involve a scmiclassical 
propagator and semiclassical vertices^3- the former is the inverse of the Hamiltonian of a Poschl- Teller problem 
in one dimension; the latter include a background dependent cubic vertex, in addition to the quartic vertex of the 
double- well. 

In this paper, we use a closed analytic form for the semiclassical propagator that has been recently obtained|l6j. 
and show that it amounts to summing up the spectral representation for the inverse of the Poschl- Teller Hamiltonian. 
Indeed, from our expression, we recover the eigenvalues and eigenfunctions of that problem. We then use that compact 
expression to rederive previous results, to extend them to lower dimensions, and to compute two-point correlations 
up to one-loop. 

The two-point vertices that we compute correspond to the self-energy at large relative distances (zero relative 
momenta). They depend on the position of those points relative to the interface, a consequence of the breaking of 
translational invariance. That self-energy is the sum of a mass (inverse correlation length) squared and a potential 
which reflects the influence of the interface on the fluctuation modes. We comment on how to compare our predictions 
for those quantities to experiments and simulations. 

The article is organized as follows: in Section II, we introduce the field theory model, and outline the derivation of 
the mean-field interface solution and of the semiclassical propagator, as well as the semiclassical expansions for the 
various generating functionals; in Section III, we compute the surface tension from the vacuum bubbles; in Section 
IV, we extract the interface profile from the graphs for the one-point Green function; in Section V, we obtain the 
self-energy from the graphs for the two-point vertex function; Section VI presents conclusions and suggestions for 
further work. Appendix A presents a detailed derivation of the semiclassical propagator, and explores its properties. 
Appendix B describes the renormalization procedure we have adopted. 



II. THE FIELD THEORY MODEL 



We consider the generating (partition) functional for a self-interacting scalar field theory model 

Z[j] = j>[Vv] cxp I - I + 1 ^ d'^xj , (1) 
whose action functional in d spatial dimensions is given by 

sM = J d'^x [\{Vipr + ^{^' - ^ir] , (2) 

where A is a dimensional coupling, ify is the vacuum value of the field (p, and j{x) is an external current. 

The model may be used to describe an interface and its fluctuations: for instance, we may associate the scalar 
field to the difference in concentration of the phases (components) of a binary system, as in a binary liquid mixture. 
Similarly, we may relate it to an Ising-like spin. The phases of the liquid mixture are described by the two degenerate 
vacua ±ipy of the model. Analogously, the two homogeneous configurations where the spins are either all up, or all 
down, may likewise be described by those two vacua. The mean-field interface will emerge as a classical solution of the 
equation of motion of the model. Its fluctuations will be captured by a semiclassical expansion around that solution. 

The classical solution of interest is the well-known kink profile, which depends on only one longitudinal coordinate 

z, 

(p{z) — ±ipyta,iih9{z) , (3a) 



= — {z~z), (3b) 

where M = (fv/\/3. The solution breaks translational invariance along the longitudinal direction as it depends on 
a position z, the point where it vanishes, which can be identified with the position of the kink. For two and three 
dimensions, the dependence of z on the transverse coordinate(s) characterizes the interface. 

The semiclassical expansion includes fluctuations of the interface in a systematic wav| l^ Il7| . Setting 



(fix) = 0{z) + X^^'^rjix) , 



(4) 
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we expand the action functional around the classical solution, regarded as a mean-field profile to be modified by the 
fluctuations, and perform a saddle-point integration to obtain the generating functional of Eq.(^) in the form of an 
infinite (semiclassical) series. 

The quadratic term in the functional expansion of the action defines the semiclassical propagator G{x,x') around 
the kink background 

[ - + _ M^)] G{x, x') - 5\x ~ x') . (5) 
Using the expression for ip, and Fourier transforming the transverse coordinates, leads to 

[-dl + P + - ^M^sech^e] G{k; z,z')^6{z- z') , (6) 

where k is the transverse momentum, and G is a hybrid momentum-position propagator. 

G can be viewed as the inverse of the Hamiltonian for a Schrodinger problem in one dimension with a Poschl- Teller 
potential U{e) = P + M"^ - {?,/2)M\ec\i^0. We may obtain an expression for G from two linearly independent 
solutions of the homogeneous version of Eq.®. It is a hypergeometric equation, but it so happens that its two 
independent solutions are hypergeometric series that terminate, as shown in Appendix A. Therefore, we end up 
with an expression for G which can be written in terms of elementary functions. Using the dimensionless quantities 
K = 2k/ M, w = (1 - tanh6')/2, and b = V4 + 

^ = ^ t-^^^' ")'^(-^' "')®("' - + ")'^(^' "')®(" ~ "')] I ' (7^) 
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In Appendix A, we show that the above form is equivalent to the usual spectral sum over eigenmodes. In fact, we use 
it to derive the eigenvalues and eigenfunctions of the Schrodinger problem. 

The cubic and quartic terms in the functional expansion of the action define the vertices of the semiclassical series 
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The cubic vertex involves the kink background. Saddle-point integration implies expanding the exponential of the 
cubic and quartic parts of the action in a power series, and then performing the resulting functional integrals, which 
reduce to powers of the fiuctuation times a Gaussian, leading to a series for the generating functional. 

The translational invariance of the model requires the introduction of a collective coordinate. The background 
solution explicitly breaks the invariance along the longitudinal direction (it is centered on z). However, as there are 
solutions for any value of z, translational invariance must be restored in the calculation by adding over all possible 
values of z. Indeed, translational invariance manifests itself through the appearance of a zero-energy eigenmode of 
the fiuctuation kernel of Eq.®. Physically, it costs zero energy to infinitesimally translate a classical solution. The 
restoration of the symmetry is accomplished by means of the Faddeev-Popov procedure, and yields a Jacobian. We 
trade integration over the zero-mode subspace for integration over the collective coordinate z. 

The generating functional, after the introduction of the collective coordinate, becomes 
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We have introduced a longitudinal infrared cutoff L, the longitudinal size of the system; the (S/XY^'^ factor comes 
from the Jacobian, while the (27r)~^/^ comes from the functional measure; A' is the determinant of [G"]"""^, with G' 
denoting the semiclassical propagator with the zero-mode subspace excluded; z'[j] = exp{i J j{x)G'{x, y)j{y)d'^xd'^y}; 
and V [if, rfi is defined by the power series expansion of the exponential in 
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SS = J {V3 + V4)d'^x. (10b) 

W[(p] is the derivative of the double-well potential U = {ip'^ — ipiy/A\ with respect to (f, computed at the kink 
solution; V3 and V4 are the vertices of Eq. (|SJ) . The Faddeev-Popov Jacobian expanded around (f yields the combination 
(S/XY/'^'P appearing in Eq.Q. 

Connected correlation functions can be derived from the free energy functional F[j] = — limi_K5o{log by 
functional differentiation with repect to j{x). A Legendre transform leads to the effective action (Gibbs) functional 
T[4>\ = F[j] + J j (j)/VXd''x, with (/>(x) = {^pix)) denoting the expectation value of the field if. Its functional deriva- 
tives with respect to (j){x) lead to one-particle irreducible {IP I) vertex functions. Connected correlations and IPI 
vertex functions can all be expressed in the usual language of Feynman diagrams, as in the explicit examples of the 
forthcoming sections. We note that a semiclassical expansion around (p (ipv) leads naturally to correlations and 
vertex functions in the interface (vacuum) sector, as we take functional derivatives at j — 0, for the former, and at 
4> = (p {'Pv), for the latter. 

Finally, renormalization is required to connect the various correlation functions to physical parameters. As we are 
interested in comparing correlations in the presence of an interface with those of the homogeneous (vacuum) phase 
(i.e., with bulk values), we will use counterterms computed in the vacuum sector. The diagrams in the kink and 
in the vacuum sector will then be related to the bulk physical parameters by means of a standard renormalization 
procedure, which is outlined in Appendix B. 



III. THE SURFACE TENSION 



The surface tension for the interface may be obtained from the difference between the free energy functionals at 
zero external current in the interface and vacuum sectors 

..lim^.lim«-«, (11) 

where A denotes the "area" spanned by the transverse direction(s). Its leading (zero-loop) value is given by the 
difference of the classical actions in the two sectors ctq = IvoiA^ooiS / \A) — 2M^/X, which is independent of the space 
dimension. 

Up to one-loop order, the free energy and effective action functionals coincide 18i], so we have 

Af[o] = r[^]-r[^.] = | + liog(^|^^) , (12) 

where A^, is the (free) determinant of the quadratic fluctuation kernel in the vacuum sector. Neglecting terms 
that vanish as {log A/ A), we have 



a = ao+ai = ao + ^lim (i^log|-) . (13) 

The (Ti contribution can be computed from the semiclassical propagator of Appendix A, as shown in reference |16| . 
Indeed, 



1 f d'^-^k , A( 



^^'^ (14) 



(fc) ' 

where A(A:) and Ai,(fc) are the determinants of [G'(fc)]^^ and [Gi,(fc)]^^, and we sum over transverse momenta. 
Exclusion of the zero-mode for fc = is guaranteed by writing 

1 1 2.,x 1 A'(fc) , A„(fc) , A'(0) 

The second term on the rhs is expressible in terms of the semiclassical propagator 

A'(k) f^^ f°° —I 
log ^7^= X '^^ J dzG{^s-z,z) , (16) 
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whereas the third is given by a similar expression with G replaced by . Then 



1 A'(fc) A;(fc) 



b{k) - 1 



b{k) + 1 



-21og[6(fc) + 2] + log 48 . (17) 



The fourth term may be computed following reference and yields log[A'(0)/At,(0)] = — log 48. Using H14() and 
p7|l , we derive the unrenormalized expression 

_ 1 f d'^-^k r [b{k) - 1] [b{k) - 2] \ 
"^-27 (2^ [fe(fc) + 1] [Kfc) + 2] / ' ^'^^ 

which is valid for fc > 0. 

In rf = 2 and d = 3, expression (|18|) needs to be regularized and renormalized. We use a cutoff in transverse 
momentum space to regularize: in d = 2, we integrate over the interval [—A, A], whereas, in c? = 3, we integrate over 
a disk of radius A. Neglecting terms that vanish as 1/A, we obtain 

3M , /2A\ / 1 3 \ , „ , 
= log — + , d = 2, (19a) 



27r \M J V4\/3 27r 
^r = -^+r^-^)M^ ^ = 3. (19b) 



47r V87r 327r 

The renormalization procedure is described in Appendix B. We adopt renormalization conditions at zero momenta 

r^'^(O) = ; Vf{Q) = , (20) 

which specify the Fourier transformed one and two-point IPI vertex functions. They lead to the renormalized result 



-i.= (^-|) , d = 2, (21a) 

^iR ^ ^ {logs - A) , d = 3. (21b) 

For d = 2, the result coincides with the correction to the kink mass obtained previously [T9l|. The results for d = 3 are 
shown, in Appendix B, to coincide with those in the literature 2, 7. 10, 12, 14]. 

We could go beyond one-loop in our expansion by including contributions which are first-order in the cubic and 
quartic semiclassical vertices. Using equations @ and (|10|) to compute F[0] to that order, we obtain the two- loop 
Fcynman diagrams and the Jacobian contribution depicted in Figure^ Such diagrams were computed in reference [iSj , 
using the spectral sum representation for the semiclassical propagator. We believe that our (resummed) expression 
for the propagator will confirm their results, and simplify the calculation, but we shall postpone that verification 
for a future publication, and concentrate on one-loop correlations in the present article. As the two-loop results 
have been compared to those obtained from Monte Carlo simulations |20| |. it is important to have an independent 
check. Presumably, the calculation will be more direct if one makes use of the compact expression for the resummed 
propagator. 

^-o ^ o o-o 

(a) (b) (c) (d) 

FIG. 1: Diagrams which contribute to the surface tension at two-loop order. 



IV. THE INTERFACE PROFILE 



The interface profile is given by the expectation value of the order parameter (j>{x) = {(p{x)), which we compute 
from the first derivative of F[j] with respect to the external current j{x), at j = 0. The calculation was carried out 
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in the kink sector using the semiclassical expansion around (p. We may express the correction to the interface profile 
in terms of the Feynman diagrams of Figure 0] As we shall show, this is equivalent to solving a one-loop corrected 
equation for (f), as was done in references It can be calculated in a much simpler and compact way using 

our semiclassical propagator. 

(a) (b) 

FIG. 2: Contributions to the interface profile. Diagram (6) represents the Jacobian term. 



Before proceeding, we return to the discussion of translational invariance. Inspection of the diagrams in Figure |21 
shows that diagram (a) is ultraviolet divergent in d = 3 if we compute it with the semiclassical propagator in the 
subspace orthogonal to the zero-mode. Indeed, the excluded zero-mode contribution behaves as 

which accounts for the divergence. The problem does not exist in lower dimensions, suggesting that d = 3 is a marginal 
dimensionQ. Were we to work with the full propagator, no ultraviolet problem would occur, but we would have an 
infrared problem for vanishing k. 

Following the interpretation of references 0,0,0] for the case of binary liquid mixtures, we take that to indicate the 
instability of a translationally invariant interface in three dimensions. To cope with this problem, we explicitly break 
translational invariance by introducing a small mass fi in the zero-mode subspace. Physically, in three-dimensional 
binary liquid mixtures that can be attributed to the action of a gravitational field, or some other pinning effect, and 
is negligible in the other subspaces (p <C M). In other physical applications, the fact that the interfaces or domain 
walls have their positions pinned down by some external effect which breaks translational invariance will be encoded 
in the dependence of /i on whatever parameter characterizes that pinning effect. 

Breaking translational invariance has the following implications for our calculation: i) we no longer need to work 
in the subspace orthogonal to the zero-mode, so that no Jacobian will emerge; ii) for the semiclassical propagator, we 
use the expression in the subspace orthogonal to the zero-mode, added to a zero-mode part with mass ^ 

Gik; = G'(fc; z, z') + f ^sMM^l , (23) 

where r^g is the normalized zero eigenmode shown in Appendix A. 

The preceding paragraph implies that the only diagram to be considered in one-loop order is the first diagram of Fig- 
ure|21 We compute it by integrating over longitudinal coordinates and transverse momenta. With the rcnormalization 
conditions defined in (|20|l . we derive the renormalized profile 

(j)R{e) = %/3A/ tanh 6* - A ( a sech^6l + /3 tanh 9 sech^g ) , d = 2 or 3 . (24) 

The profile has the functional form above, and only the coefhcients change with the spatial dimension. That is a direct 
consequence of the form of the semiclassical propagator. The computed coefficients are a — {2tt - 3V3)/(167rM), 
P = -{6^/3 + 4:TT - 3%/37rAf/^)/(247rM), for d = 2, and a = 3V3(log3 - l)/(327r), /? = (V3/167r) x log(4M/3At), for 
d = 3. It should be remarked that, for d = 3, our calculation will coincide with that of reference for a judicious 
choice of rcnormalization conditions. The latter reference confronted its findings with experimental results_l_5j, being 
compatible with the data available at the time. Our new results for d = 2, depicted in Figure O illustrate the 
dependence of the profile on the ratio /i/M. Obviously, the lower that ratio, the more striking the effect will be. 
Those results can be tested experimentally by studying the interface of thin films of binary mixtures|2lj. for instance. 
More recent applications, such as the ones involving Bose-Einstein condensates, could also be used in experimental 
checks in d = 2 and d = 3 [2^. 

As we have already remarked, it is only in d = 3 that we are forced to break translational invariance. In d = 2, 
a translationally invariant mean-field interface solution is stable, so that we may compute its fluctuations by using 
G , and including diagram (b) of Figure El However, in the limit of large transverse "area (in the present case, 
a length), that contribution is negligible. The result we obtain has the functional form presented in (|24|) with 
a = {2tt — 3-\/3)/(167rM) and P = — (6-\/3 + 47r)/(247rA:f) (which is equivalent to taking /i ^ oo in the expression for 
(3). That is to be compared with the profile obtained previously, without translational invariance. 
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FIG. 3: Profile diagrams in d = 2 for A = 0.1. The sofid line is the kink configuration. The dashed curve correspond to 
lijM = 0.01 and the dotted one to [ijM = 0.001. 



V. THE TWO-POINT CORRELATIONS 

We may take a second derivative of the free energy functional F with respect to the external current j to obtain the 
connected two-point function G^c^ . That leads to the Feynman diagrams of Figure^ Just as in the previous section, 
one has to omit the diagram coming from the Jacobian, whenever translational invariance is broken. 



(o) (6) (c) (d) 

FIG. 4: Diagrams which contribute to the two-point function up to one-loop order. 

(2) 

The calculation of the two-point function Gc up to one-loop order is rather involved, so we shall postpone it for a 
future publication. Instead, we will concentrate on mean-field results obtained from its lowest order expression, which 
is given by our semiclassical propagator, and on one-loop corrections to F*^^). 

As our classical interface profile depends on a collective coordinate z, the hybrid Fourier transform G(fc; z, z') 
depends on both z and z', not just on their difference. Introducing center-of-mass and relative coordinates 



R 
P 



z + z 
" 2 
z-z' 



(25a) 
(25b) 



we may reexpress G in terms of R and p. The expression has a particularly simple form ior k — 

^, ^ e^-^POo(M/^if - 4 - MIp) + 8e^P cosh(Mi?) + cosh{2MR) , 

G(R,p,0) = 5 B(p) -I- [p < 

M(l + eA^P + 2e*W2 cosh(M-R)) 



-P) 



(26) 



If we now perform a Fourier transform in the p-coordinate, we obtain a function G depending on R, kp and k. 
Setting kp — and fc = amounts to integrating over all relative coordinates. 
One may define a susceptibility x ^ 



1 /-^/^ 

X= lim - / dRG{R; 0,0) 



(27) 



Likewise, one may exclude the lowest mode and define x' using G' in the previous formula. For both cases, we obtain 
X = 1/M^, just as in the vacuum sector. That is a consequence of the fact that the two lowest modes are localized, 
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whereas the continuum ones behave asymptotically as plane waves of mass M . Accordingly, the exponential decay 
of our propagator as the relative distance becomes large is of the form ex'p{—M^p'^ + p^), where pr is the relative 
transverse coordinate. The correlation length is thus set by 1/M, being independent of i?, even when we give a small 
mass /i to the lowest mode. In principle, the susceptibility and the correlation length could depend on the position 
of the two points with respect to the interface, i.e., on K. However, the phases on cither side of the interface are 
degenerate in our model. They have the same correlation length 1/Af. It is then natural that M should set the scale. 
A model wherein the co-existing phases could have different masses (inverse correlation lengths) would probably lead 
to position-dependent quantities. 

The calculation of the two-point vertex F^^^ up to one-loop order involves fewer integrals than that of its inverse 
d'c'^. Setting 

r(2)(a;,a;') = [G(2)]"\x,x') = G~^{x,x') + Y.(x,x') , (28) 

where the first term on the rhs is the inverse of the semiclassical propagator, then S(x,a;') will be the contribution 
to the self-energy from the kink sector, and can be identified with two of the diagrams of Figure 0] without external 
legs, as shown in Figure EKdiagram (£) is obtained as a combination of F^^' and F^^^ with G~^). 




(a) (6) 

FIG. 5: Diagrams which contribute to the self-energy up to one-loop order. 

We should stress that the semiclassical propagator itself already gives a contribution Esc to the self-energy when 
compared with the vacuum sector 

G-\x,x') ^G;;\x,x') + i:scix,x') , (29) 

where is the inverse of the free propagator. Equation H29|l defines S^c , just as equation (|28|l defines E. Fur- 
thermore, the diagrams of Figure |S1 when computed with the free propagator, and with ip^ instead of ip at the cubic 
vertices, will yield E„(x,x'), defined by 

T<^Hx,x') = [G<^^r\x,x')^G-\x,x') + E,ix,x') . (30) 

Equations l|2HIl, ^ and ^ lead to 

AT'-^\x,x') = T'-^\x,x') -T'^Hx^x') = T.sc{x,x') + T,{x,x') - E„(x,a;') = AE , (31) 

the difference in self-energy between kink and vacuum sectors. 

As in the case of the two-point function, the hybrid Fourier transform AE(A:; z, z') will depend on both z and z', not 
just on their difference, so that we may reexpress AE in terms of R and p. Performing, as before, a Fourier transform 
in the p-coordinate, we obtain a function AE depending on R, kp and k 

AT.{R-kp,k) = f^''\R-kp,k)-T^;^\R-kp,k) . (32) 

In order to understand the physical meaning of the self-energy in the present situation, it is instructive to consider 
the lowest order term defined in Ea. H29|) . From Eq.(|SJl, we have 

G'\k; z, z')^[- dl + P+ ~ -Mhech^e]6{z - z') , (33) 

with the 9{z) previously defined. In terms of R and p, we have 

Ese(i?, p)=[- ^Mhech^9{R + ^)] Sip) . (34) 
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Integrating over the p-coordinate, i.e., taking the Fourier transform at kp = 0, yields a i?-dependent potential Vsc{R) = 
— (3/2)M^sech^6'(i?), that vanishes at infinity. Therefore, the lowest order contribution to T^^\R]kp,k) at zero 
relative momenta is given by the mass squared plus a fluctuation potential that vanishes as i? — > oo. The first-order 
contribution can be split likewise: the constant term as i? — > cx) can be interpreted graphically as coming from 
the diagrams of Figure computed at zero relative momenta with vacuum propagators, and with (p replaced by 
in diagram (6). Those are exactly the corrections to the mass (inverse bulk correlation length) squared. Our 
renormalization condition l|20(l guarantees that the constant term is just the renormalized mass squared; on the other 
hand, the contribution to the fluctuation potential, which vanishes as i? — *■ oo, is given by AI](i?; 0, 0). The fluctuation 
potential reflects the presence of the interface. 

Using our propagator to compute the diagrams, we obtain for diagram (a) of Figure |S1 

^ sech^ei ^%/3 tanh^ ^ i ~ ^) ^^^^^^ ' d = 2 , (35a) 

sech^ei ^(log 3) tanh^ e + log ^— ^ sech^6l^ , d = 3 . (35b) 

The calculations were done using Mathematica. The integrations for diagram (6) of Figure |S1 had to be performed 
numerically, so that we have computed it for several values of R in order to draw the curve of the modified potential 
illustrated in Figure El As before, we have adopted the renormalization conditions in H2U|) . The modified potential 
is very sensitive to the value of the ratio /i/M, especially in the case d = 2. A comparison between Figures El and O 
suggests that it is probably easier to measure the effects of fluctuations on the potential than on the interface profile. 

Radiation scattered through the interface may be used to probe structure factors, which ultimately measure the 
two-point function. A numerical evaluation of &^\R]kp,k) would then allow a direct comparison with data from 
binary liquid mixtures, or from more recent applications, especially those involving Bose-Einstein condensates. As 
we have said before, that calculation is feasible, but rather involved. Since r'^'(_R; fcp, k) measures the change in free 
energy (which coincides with the effective action to first-order) as the profile changes, one might hope to have a direct 
test of our computation by measuring those changes for well-separated points at different values of R. Alternatively, 
Monte Carlo computations might be used as a test. 



0,5 p 
:^ 



c 

p -1 

Oh 




FIG. 6: Potential as a function of R for fi/M = 0.1 and A — 0.1. The solid line stands for Vsc{R)', the other hnes correspond 
to the modified potential in d = 3 (dashed) and d — 2 (dotted). 



VI. CONCLUSIONS 



We have cast all previous results for the correlation functions of interface fiuctuations in the unified framework of 
a semiclassical expansion. Besides making contact with the existing literature, we have extended previous results to 
a lower dimension, making use of the closed analytic form for the semiclassical propagator. We have also computed 
susceptibilities and two-point correlations, which may eventually be checked experimentally. 

It is important to note that our technique for resumming the spectral representation for the semiclassical propagator 
might be of use in other contexts, as long as we can reduce the differential equation for the propagator to an ordinary 
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one. In the present case, an additional simplification came from the fact that the hypergeometric series involved 
terminated. In particular, spherically symmetric backgrounds would be natural candidates to be investigated. 

We should emphasize that with our basic ingredient, the semiclassical propagator, we have reduced the calculation 
of physical quantities to computing Feynman diagrams, whose propagators and vertices carry information about the 
background solution. 

Many systems of interest may profit from the semiclassical treatment that we have presented. Thus, computing 
correlations in supersymmetric models is clearly a direction for future work. Likewise, correlations for mixtures of 
Bose-Einstein condensates separated by an interface are certainly worth pursuing. For this latter example, existing 
experimental techniques|22 may open up a host of possibilities for experimental tests and checks. Also promising 
are the possibilities of comparison with experimental data in the more traditional binary liquid mixtures. We should 
emphasize that our treatment allows for a complete and separate treatment of capillary (those in the zero-mode 
subspace) and noncapillary waves. Reflectivities and form factors extracted from scattered radiation are the physical 
quantities to be measured for comparison. Monte Carlo simulations might also be used as a test. 

Finally, we may hope to calculate other correlations of interest to experimentalists, as long as their defining Feynman 
diagrams lead to tractable integrals. As mentioned before, the two-loop calculation done previously for the surface 
tension may serve as a test of the simplification introduced by our expression for the propagator. 
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APPENDIX A 

The homogeneous version of Equation ^ in the text can be written in terms of dimensionless variables 

+ sech^e] = , (Al) 

where u = (1 — tanh0)/2, as before. Defining </> = (cosh^?)^'' the function F satisfies a hypergeometric equation 

fJ^F dF 
u{l ~u) — +[{b+l)- 2(6 + l)u] — + [6 - 5(6 + 1)] , (A2) 

whose general solution is 

F{u)^cx[2Fx{h-2,b + i;l + h;u)]+C2[u-\Fx{i,-2-l~h;u)] , (A3) 
where 2-^1 (^, B] C; u) is a hypergeometric function. The identity 

2Fi(A, B; C- u) = {l- uf-^-^ 2Fi{C - A,C - B;C;u) , (A4) 

yields 

F(m) = ci(l-u)-^ [2Fi(3,-2;l-f 6;m)] +C2U-^[2Fi(3,-2;l-6;?x)] . (A5) 

Both series terminate. The two solutions in the linear combination of ljA5|l correspond to (f){b, u) and (j>{—b, u) of 
equation lj7b|l of the text. 

From the solutions of the homogeneous equation, one constructs the semiclassical propagator by a standard 
procedure |l^, which leads to H7al) . That expression can be rewritten as 



G = f{b, u) f{-b, u') e{z - z') + f{-b, u) f{b, u') Q{z' - z) . (A6) 

The latter expression appears as one of the terms of the integral 

dq 4'{q,z)r{q.z') 
J-oo^T^ q^ + k^ + AP ' 
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where ^p{q,z) = e"'^/(— Indeed, its Q{z — z') part has poles in the upper-half of the complex g-plane at 
iyk^+M^ , ~iM and +iM , with residues given by 

^_Vfe^+Af2(._,') ^^^^ u)f{-b, u')/^^i\Jp + M^ , (A8a) 

- 12M m(1 - u)u'(l - u') /A-KiP , (A8b) 



3M 



-f2 



6Mv/u(l-u)(l - 2u) Vu'(l-u')(l - 2M')/47ri(fc2 + ) , (A8c) 



respectively. The Q{z — z') part has similar contributions. Combining them, we derive 

2 y fc2 1^ 2 y fc2 + 3M 

which leads to the spectral representation for G 

^ AI f f]Q{z)f]Q{z') fii{z)fji{z') \ dq i}{q,z)i:*{q,z') 

2 1 F fc2 + 3Mi / 2tt ^2 + p + ' ^ ^ 

where one clearly identifies the eigenvalues and the eigenfunctions for the Poschl- Teller problem at hand. Using the 
variable u, we have 

Ao - ; fjoiu) = 2 V3uil-u ) , 

Ai=^2+3Mi. fj,{u) = ^6u{l~u){l-2u) , (All) 

A, = q2 + fc2 ^ ^.f2 . f^^^u)^e"'' f{-2iq/M,u) . 

Having shown that our semiclassical propagator does correspond to the usual spectral representation, we now 
investigate its limit when k — > (6 — > 2) in the subspace orthogonal to the zero eigenmode. That limit is required for 
the calculations of Section V. Unfortunately, the expression for the limit which appeared in reference |0| is wrong (it 
is not orthogonal to the zero- mode subspace) . We profit from this occasion to exhibit the correct expression 



w(l-u') 



u'{l - u) 



(A12) 

-^u{l-u)u'{l-u')le{u' -u) + {u<-^u'), 

where 

Giu, u') = (1 - u)^ + 6u'(l - u)2 -I- 6u'^(l - m) + u'^ . (A13) 

The corrected expression for G is indeed orthogonal to the zero mode subspace, as can be verified in a straightforward 
calculation. 

APPENDIX B 

In this Appendix, we shall outline the renormalization procedure adopted in the text. We start from the effective 
action functional, up to one-loop order 



Aicj^ix)] ^ \mx)] = smx)] + ^ log (^^) , (Bi) 



written in terms of renormalized parameters, and add to it counterterms, in order to obtain a renormalized expression 

(B2) 
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In the formulae above, (j>{x) = {^p{x)) is the expectation value of the field. The renormalization constants Ci, C2, C^l 
and C3T are associated to mass, coupling, and longitudinal and transverse wave function renormalization. They will 
be fixed by renormalization conditions in the vacuum sector. 

Functional derivatives of l|B2l) with respect to (f>{x) lead to the n-point vertex functions. Derivatives taken at 
0(a;) = 'fix) yield vertices in the kink sector. As ip satisfies the equation of motion, for the renormalized one-point 
function in the kink sector, we obtain 



r 



(i)r 



A 



(B3) 



where all quantities are taken at a given point x, so that G = G{x, x) is the semiclassical propagator at coincident 
points. For the two-point vertex, we have 



A 



A. 



T)j'[ip;xi,X2] =G '^(xi,X2) + - 2:2) G - -<^(a:i)G(a;i, a:2)^(2:2)G(a;2, a:i) 

~ [Gi + G2(3^2 „ _ c^^q2 _ C3tV|] S^^Hxi - X2) , 



(B4) 



where again (p — (p{xi), and G — G{xi,xi). 

Expressions for the vacuum sector are obtained from functional derivatives at ify. In the formulae above, that 
amounts to replacing ip and G with ip^ and G«, respectively. One may derive expressions for F^' and F^\ as well. 
The translational invariance in the vacuum sector makes it convenient to go to momentum space. Furthermore, we 
shall adopt zero-momentum renormalization conditions. If we define the Fourier transformed vertices as 



(B5) 



we arrive at the following relations at zero-momenta 



-(fy I Gv{k) - Cl(fv , 



^R\p^ = 0) = Af2 + I Gy{k) ^ I Gy{k) - Gi - 2G2(^: 



\d[pl 



= 1 - 



2 d[pi] 
2 d[pl' 

3A 



G,(fc)G„(fc + p) 
G.(fc)G„(fc+p) 



G 



3L 



Pi=0 



Pi=0 



G3T , 



^R\p^ = 0) = \-f I Giik) + 6\^i I Giik) - 3\ft I Glik) ~ 6G2 



(B6a) 
(B6b) 

(B6c) 

(B6d) 
(B6e) 



where all integrals are calculated with a cutoff in transverse momentum space. If we define 



/„KA) 



d-^-^kr r°° dk, 



27r 



G^ik) 



(kx and fc^ amount to transverse and longitudinal momentum, respectively) relations (|B6|I become 



(B7) 



r^'^(p. - 0) = + A) - ^/2(rf, A) - Gi - 2G2^^ 



dV 



dm 
d[pl] 



= 1 - 



= 1 - 



ih{d, A)-AM'^h{d, A) 



d+2 



h{d + 2,A) 



G 



3L 



hid, A) 



~dh 



^R\p^ = 0) = A - y/2(d. A) + 6A(^^/3(d, A) - 3A^::/4(d, A) ~ 6G2 , 



d+2 
2 ; 



/4(rf+2,A) 



-G 



3T 



(B8a) 
(B8b) 

(B8c) 

(BSd) 
(BSe) 



13 



where is the usual d-dimensional soHd angle. 

The renormalization conditions (|20|l for the vacuum sector which were adopted in the text lead to the determination 
of the constants: Ci = A/i (d, A)/2, C2 — —A/2 {d, A) /4, C3L = C3T = 0. Using such values in the A-regulated Fourier 
transformed expressions ljB3|) and l|B4l) at zero momenta cancels the ultraviolet divergences as A ^ c». 

As a consistency check on our procedure, we have used the renormalization conditions 



that were adopted in reference in order to calculate corrections to the kink mass. The results we have found 
coincide with those of 0, illustrating that the calculation of the determinant via the semiclassical propagator is not 
afflicted with mode-counting ambiguities. That allows us to use a simple momentum cutoff regularization, which is 
more directly related to condensed matter phenomenology. 



[1] T. W. B. Kibble, in Topological Defects and the Nonequilibrium Dynamics of Symmetry Breaking Phase Transitions, Les 
Houches, France (1999). 

[2] A. Rebhan and P. van Nieuwenhuizen, Nucl. Phys. B 508, 449 (1997); H. Nastase, M. A. Stephanov, P. van Nieuwenhuizen 
and A. Rebhan, Nucl. Phys. B 542, 471 (1999); A. Rebhan, P. van Nieuwenhuizen and R. Wimmer, New J. Phys. 4, 31 
(2002). 

[3] T. Bhattacharya, A. Gocksch, C. Korthals Altes and R. D. Pisarski, Phys. Rev. Lett. 66, 998 (1991); Nucl. Phys. B 383, 
497 (1992). 

[4] D. T. Son and M. A. Stephanov, Phys. Rev. A 65, 063621 (2002). 
[5] A. Parnachev and L. G. Yaffe, Phys. Rev. D 62, 105034 (2000). 

[6] E. Farhi, N. Graham, P. Haagensen and R. L. Jaffe, Phys. Lett. B 427, 334 (1998); N. Graham and R. L. Jaffe, Nucl. 
Phys. B 544, 432 (1999); tbid. 549, 516 (1999); E. Farhi, N. Graham, R. L. Jaffe and H. Weigel, ihid. 585, 443 (2000); 
ihid. 595, 536 (2001); N. Graham, R. L. Jaffe, M. Quandt and H. Weigel, Phys. Rev. Lett. 87, 131601 (2001); N. Graham, 
Phys. Lett. B 529, 178 (2002). 

[7] D. Jasnow, in Phase Transitions and Critical Phenomena, vol. 10, edited by C. Domb and M. S. Green (Academic Press, 
New York, 1986), p. 270. 

[8] D. J. Wallace, in Recent Advances in Field Theory and Statistical Mechanics, Les Houches, 1982, edited by J.B. Zuber and 
R. Stora (North Holland, 1984). 

[9] J. Rudnick and D. Jasnow, Phys. Rev. B 17, 1351 (1978). 
[10] G. Miinster, Nucl. Phys. B 324, 630 (1989). 
[11] J. Rudnick and D. Jasnow, Phys. Rev. Lett. 41, 698 (1978). 

[12] G. Miinster, Nucl. Phys. B 340, 559 (1990); S. Klessinger and G. Munster, Nucl. Phys. B 386, 701 (1992). 

[13] P. Hoppe and G. Miinster, Phys. Lett. A 238, 265 (1998). 

[14] E. Brezin and S. Feng, Phys. Rev. B 29, 472 (1974). 

[15] E. S. Wu and W. W. Webb, Phys. Rev. A 8, 2065 (1973). 

[16] C. A. A. de Carvalho, Phys. Rev. D, 65, 065021 (2002). 

[17] C. A. A. de Carvalho, R.M. Cavalcanti, E.S. Fraga and S.E. Joras, Annals Phys. 273, 146 (1999). 

[18] D.J. Amit, Field Theory, The Renormalization Group and Critical Phenomena (World Scientific, Singapore, 1984). 

[19] R. Rajaraman, Phys. Rep. 21C, 227 (1975); Solitons and Instantons: an Introduction to Solitons and Instantons in 

Quantum Field Theory (North-Holland, Amsterdam, 1987-89), and references therein. 
[20] M. Hasenbusch and K. Finn, Physica A 245, 366 (1997). 
[21] K. Binder, Phys. Rev. Lett. 77, 1318 (1996); J. Stat. Phys. 95, 1045 (1999). 
[22] F. Dalfovo, S. Giorgini, L. P. Pitaevskii, and S. Stringari, Rev. Mod. Phys. 71, 463 (1999). 




(B9) 



